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The generalized master equation with two times, introduced in [1,2], applies to 
the problem of diffusion in an time-dependent (in general inhomogeneous) external 
field. We consider the case of the quasi Fokker-Planck approximation, when the 
probability transition function for diffusion (PTD-function) does not possess a long 
tail in coordinate space and can be expanded as the function of instantaneous dis- 
placements. The more complicated case of the long tails in PTD will be discussed 
separately. 

I. INTRODUCTION 

The models of continuous time random walks (CTRW) [3], when some objects can jump 
from one point to another in inhomogeneous (in general) media and stay during some time 
in these points before the next usually stochastic jump, are important for the solution of 
many physical, chemical and biological problems. Recently these models have been applied 
also in economics and social sciences (see, e.g., [4-6]). Usually the stochastic motion of 
the particles leads to a second moment of the density distribution that is linear in time 
< r 2 (t) >~ t. Such type of the diffusion processes play a crucial role in plasmas, including 
dusty plasma [7], nuclear physics [8], neutral systems in various phases [9] and in many 
other problems. However, in many systems the deviation from the linear in time dependence 
of the mean square displacement have been experimentally observed, in particular, under 
essentially non-equilibrium conditions or for some disordered systems. The average square 
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separation of a pair of particles passively moving in a turbulent flow grows, according to 
Richardson's law, with the third power of time [10]. For diffusion typical for glasses and 
related complex systems [11] the observed time dependence is slower than linear. These two 
types of anomalous diffusion obviously are characterized as superdiffusion and subdiffusion. 

The generalized master equation for the density evolution, which describes the various 
cases of normal and anomalous diffusion has been formulated in [1,2] by introduction of the 
specific kernel function (PTD) depending on two times W(r, r', r, t — r), which connects in 
a linear way the density distributions of the stochastic objects (or particles) / for the points 
r' at moment r and r at moment t. The approach suggested in [1,2] clearly demonstrates 
the relation between the integral approach and the fractional differentiation method [12] 
and permits to extend (in comparison with the fractional differentiation method) the class 
of sub- and superdiffusion processes, which can be successfuly described. On this basis in 
[2] the different examples of superdiffusive and subdiffusive processes were considered for 
the various kernels W and the mean-squared displacements have been calculated. Recently 
the idea of the generalized master equation with two times [1,2] for diffusion in coordinate 
space has been recently used in [13] for the calculation of average displacements in the case 
of a time- dependent homogeneous external field. 

This paper is motivated by the necessity to describe in more detail the influence of time- 
and space dependent external fields on the continuous-time random walks. The equation 
formulated in [1,2] is appropriate for this purpose and offers the opportunity for consideration 
of CTRW for both cases: long-tail space behavior of the PTD function, as well as for the 
fast decay of PTD function in coordinate space, when the Fokker-Planck type expansion is 
applicable. For simplicity in this paper we consider only the last case. 
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Equation can be represented in an equivalent form, more similar to the structure of the 
Fokker-Planck equation, where the initial condition is absent: 



II. GENERALIZED MASTER EQUATION 



Let us start from the generalized master equation with two times [1,2]: 




df(v,t) _ d 



dr / dr' {W(r, r', r, t - r)f[r', r) - W(v', r, r, t - r)/(r, r)} . (2) 



dt dt 
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or 

df(r,t) 

at 

where the function P(r, r', r,t — r) is given by: 



dr J dr' {P(r, r', r, t - r)f(r', r) - P(r', r, r, t - r)/(r, r)} , 



(3) 



P(r, r', r, t - r) = 2W(r', r, r, t - r)5(t - r) + J^(r', r,r,t- r) (4) 

The argument t— r describes the retardation effects, which can be connected in the particular 
case of multiplicative PTD function W(r, r', r,t — r) = W(r, r', r)x(t — t) with, for example, 
the probability for the particles to stay during some time at the fixed coordinate before 
transition to the next point. An equation with retardation, where W function depends only 
on one time argument t — r has been suggested first in [14] and applied in [15] to the case of 
the multiplicative representation of the PTD function. In general W is not a multiplicative 
function in the sense mentioned above and, what is more important, is a function of two times 
t and t — r [1]. It is necessary to mention that the closed form of the equation for the density 
distribution is an approximation. In some cases the exact solution for density distribution 
can be found (see e.g. [16], [17]), when the closed equation for density distribution does not 
exist or gives an approximate result. Nevertheless, in many practical situations Eq. ((TJ or 
(Jl| are sufficiently exact and permit to describe various experimental data. 

Let us consider the nature of appearance of the two time arguments in the generalized 
master equation Eq. in the case of a time-dependent external force F(r, t). To simplify 
the consideration we can investigate the case of fast decay of the kernel W(r, r', r,t — r) = 
W(u, r, T,t — r) as a function of u = r — r', when an expansion in the spirit of Fokker-Planck 
can be applied. In this case Eq. ((TJ takes the form [1,2]: 



/(r,f) = /(r,f = 0) + fdr-^r 



d 

A a (r, r, t - r)/(r, r) + — (B aP (r, r, t - r)f(r, r) 



,(5) 



where the functions A a (r, r, t — r) and B a p(r, r,t — r)/ ff (r, r) are the functionals of the PTD 
function (the indexes are equal a,/3 = x s in s-dimensional coordinate space): 



A a (r, r, t — r) = J d s uu a W(u,r,r,t - r) (6) 

and 

B a p(r,r,t - t) = i J d s uu a upW(vL,r,T,t-T). (7) 



4 



Eq. (jSJ) can be rewritten naturally in the form, corresponding to Eq. (J2J), but now for the 
Fokker-Planck type approximation: 



df(r,t) d [ l d 



d 

A a (r, r, t - r)f(r, r) + — (B a(3 (r, r, t - r)f(r, r)) 



dt dt J dr a 

We suggest, that the PTD function is independent of f(r, t), therefore the problem is linear. 

III. INFLUENCE OF THE EXTERNAL FIELDS 

One of the main sources of inhomogeneity is an external field, which also provides the 
prescribed dependence of the PTD function on r. Other wards we can suggest, in the 
considered particular case, that the dependence of VF(u, r, r, t — r) on the arguments r, r is 
connected with a functional dependence on an external field: 

W(u, r, r, t - t) = W(u, t - r; F(r, r)). (9) 

If the external fields are absent the PTD function is a function of the modulus u = u, which 
means that A a = and B = SapB^it — r) with: 

B (t - r) = I ' d s u u 2 W (u, t-r). (10) 

ZS J 

For relatively weak external fields the functional (JHl) can be linearized in the external field 

W(u, t - r- F(r, r)) = W (u, t-r)+W x {u,t- r)(u ■ F(r, r)). (11) 

The functions Wo(u, t— r) and Wi(u,t — r) are equal to W(u,t—r; F = 0) and the functional 
derivative 5W(u,t — r; F(r, r))/5(u • F(r, r))| F=0 respectively. Then the functions A a and 
B a p take the form 

A a (r, T,t-r) = ~F Q (r, r) J d'wfW^u, t - r) = F a (r, r)L(t - r), (12) 
where L(t — r) is given by 

L{t-r) = - s J d s uu 2 W!{u, t-r). (13) 

and 

Bap(T,T,t-T)=5 a pB (t-T). (14) 
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The generalized diffusion equation Eq. (jHJ) takes the form 

= j t j^dr [L(t - r) V(F(r, r)/(r, r)) + B (t - r) A/(r, r)] , (15) 

In general this equation contains two different functions -Bo and L depending on the argument 
t — t. If the functional W^u, t — r; F(r, r)) is multiplicative, namely, W^u, t — r; F(r, r)) = 
iy(u; F(r, r))x(t — r) Eq. (|15p can be simplified: 

df(r,t) 



j t J* drxit - r) [DAf(v, r) - 6V(F(r, r)/(r, r))] , 



dt 

Here 6 and D are the constants, determined by the relations: 



(16) 



b = -~y rfWiy^M) (17) 

with W^(u) = <W(u; F(r, r))/5(u ■ F(r, r))| F=0 and 

D = — [ d s uu 2 W (u). (18) 

The physical sense of the multiplicative structure of the functional W is that independence of 
the time delay of the random walkers is independent of the external field. The dimensionless 
function x{t) in this simple case is connected with the hopping-distribution function ip(t) = 
\ip*(\t) introduced in the master equation by Scher and Montroll [15]. The value A = 1/Vq 
is the characteristic waiting time for the hopping-distribution. Laplace transformations of 
these functions x( z ) an d V ; *(- 2 ) relate them as follows 

y(z) = r{z) (19) 

For the exponential hopping-time distribution ip(t) = \exp(—\t), where A = 1/tq (tq is the 
characteristic waiting time) we have ip*(z) = 1/(1 + z), x( z ) = 1/z and x(t) = xfa) = 1- 
In this case Eq. f|16j) is reduced to the usual diffusion equation in an external field with the 
diffusion coefficient D and mobility b: 

= DA f(r, t) - bV (F(r, t)f(v, t)) . (20) 
IV. CONCLUSIONS 



We show that the generalized master equation with two times, which have been introduced 
in [1,2] can describe the influence of inhomogeneous and time-dependent external fields on the 
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diffusion processes. Linearization of the general master equation in the external field leads 
to essential simplifications. In this case the diffusion processes depend, in general, on two 
different functions of time, which describe retardation due to the finite time of occupation 
and transferring particles in space in the presence of the external field. Relations with 
simpler models are established. The consideration in the present paper give the opportunity 
to consider a wide class of the problems of normal and anomalous transport in external fields 
on the basis of generalized master equations with two times. 
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